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ABSTRACT 

An  extensive  quantity  of  airload  measurements  was  ob¬ 
tained  for  a  pressure-instnimented  model  of  the  BO- 
105  main  rotor  for  a  large  number  of  higher-harmonic 
control  (HHC)  settings  at  Duits-Nederlandse  Wind 
Tunnel  (DNW).  The  wake  geometry,  vortex  strength, 
and  vortex  core  size  were  also  measured  through  a 
laser  light  sheet  technique  and  LDV.  These  results  are 
used  to  verify  the  BVI  airload  prediction  methodolo¬ 
gies  developed  by  AFDD,  DLR,  NASA  Langley,  and 
ONERA.  The  comparisons  show  that  an  accurate  pre¬ 
diction  of  the  blade  motion  and  the  wake  geometry  is 
the  most  important  aspect  of  the  BVI  airload  predic¬ 
tions. 

1  INTRODUCTION 

Although  the  current  civil  helicopter  has  been  used  ex¬ 
tensively  because  of  its  excellent  hover  and  low  speed 
forward  flight  capability,  the  civil  helicopter  fleet  con¬ 
tinues  to  be  too  noisy  for  community  acceptance,  pre¬ 
venting  it  from  reaching  its  full  potential  of  utilization. 
One  of  the  major  noise  sources  comes  from  the  rotor 
blade  cutting  through  its  own  wake.  This  phenomenon 
is  known  as  blade- vortex  interaction  (BVI).  It  occurs 
mostly  during  descent  flight  for  landing. 

Over  the  last  decade,  many  experimental  studies  from 
full-scale  flight  tests  (Refe.  1,  2)  to  model-scale  tests 
(Refs,  3,  4,  5)  have  been  performed  to  investigate  this 
BVI  phenomenon.  It  was  foimd  that  the  BVI  noise 
can  be  scaled  (Refs.  6,  7),  and  expensive  flight  tests 
can  be  simulated  by  model  rotor  tests.  The  develop¬ 
ment  of  BVI  prediction  methodology  lagged  behind 
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the  experimental  studies  since  no  comprehensive  ro¬ 
tor  code  had  enough  azimuthal  (or  time)  resolution 
to  predict  the  BVI  event.  However,  several  analytical 
methods  (Refe.  8,  9)  were  developed  with  limited  suc¬ 
cess  as  compared  with  measured  blade  surface  pressure 
distributions  and  acoustic  signatures.  Refe.  8  and  9 
indicate  that  the  wake  geometry  and  vortex  strength 
are  two  major  parameters  affecting  the  BVI  noise  sig¬ 
natures.  The  information  of  the  wake  geometry  and 
vortex  strength  were  not  measured  in  most  tests. 

An  international  cooperative  test,  called  the  Higher- 
harmonic-control  Aeroacoustic  Rotor  Test  (HART), 
was  conducted  by  AFDD,  DLR,  DNW  and  ONERA 
at  Duits-Nederlandse  Wind  Tunnel  (DNW)  to  mea¬ 
sure  not  only  the  blade  surface  pressure  distributions, 
the  blade  deformation,  and  the  acoustic  signatures, 
but  also  the  wake  geometry  using  Laser  Light  Sheet 
(LLS)  technique  (Ref.  10)  and  the  vortex  strength  us¬ 
ing  Laser  Doppler  Velocimetry  (LDV)  method  (Refe. 
11, 12).  The  measured  data  will  be  used  to  validate  the 
prediction  methodology  and  to  set  the  future  direction 
for  code  improvement. 

Analytical  predictions  of  BVI  events  by  AFDD,  DLR, 
NASA,  and  ONERA  were  performed  before  the  HART 
test.  The  predicted  wake  geometry  served  as  a  guide 
for  the  LLS  and  LDV  teams  to  locate  the  desired  vor¬ 
tex.  Thus,  the  wind  tunnel  simulations  of  a  6®  descent 
flight  case  without  HHC  and  two  similar  cases  with 
HHC  are  selected  for  comparison.  One  is  denoted  as 
the  baseline  case.  The  second  one  is  the  low  noise 
case  and  the  third  one  is  the  low  vibration  case.  The 
results  for  comparison  include  the  blade  flapping  mo¬ 
tion,  elastic  deformation,  blade  surface  pressure  distri¬ 
bution,  sectional  lift,  and  wake  geometry  at  selected 
azimuthal  angles. 
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2  PREDICTION  METHODOLOGY 

Before  the  HART  test,  a  prediction  team,  consisting 
of  researchers  from  US  AFDD,  German  DLR,  NASA 
Langley,  and  French  ONERA  was  formed.  Each  mem¬ 
ber  of  the  prediction  team  predicted  several  cases  to 
assist  the  test  team  in  finding  the  BVI  locations.  The 
prediction  methodology  developed  by  each  organiza¬ 
tion  has  been  discussed  in  detail  (Ref.  9).  A  brief 
description  of  each  prediction  method  is  given  in  the 
following. 

To  predict  the  blade  loading,  AFDD  uses  a  comprehen¬ 
sive  analytical  model  of  rotor  aerodynamics  and  dy¬ 
namics  (C  AMR  AD/ J  A)  to  determine  the  vortex  tra¬ 
jectory,  vortex  strength  and  blade  partial  inflow  an¬ 
gles.  This  information  is  then  used  as  an  input  to 
FPR,  a  3-D,  unsteady  full  potential  code,  to  deter¬ 
mine  the  aerodynamic  surface  pressures.  Similarly, 
DLR  uses  the  S4  code  which  consists  of  three  mod¬ 
ules;  namely,  an  unsteady  aerodynamic  model,  a  wake 
model,  and  an  aeroelastic  model.  These  three  mod¬ 
els  are  controlled  by  the  rotor  trim  algorithm.  The 
aerodynamic  model  is  a  lifting  line  method  with  the 
Leiss  dynamic  stall  algorithm.  The  wake  model  uses 
a  prescribed  and  distorted  wake  developed  by  T.  Bed- 
does  (Ref.  13).  Uncoupled  flap,  lead-lag,  and  tor¬ 
sional  modes  are  used.  At  NASA  Langley,  a  compre¬ 
hensive  rotor  code  CAMRAD.MODl  and  HIRES  or 
a  full-potential  CFD  code,  FPRBVI  are  used  for  the 
blade  load  prediction.  The  numerical  methods  devel¬ 
oped  at  ONERA  combine  R85  and  METAR  codes  to 
get  the  trim  solution,  MESIR/MENTHE  to  calculate 
the  wake  geometry  and  ARHIS  to  obtain  the  blade 
loading. 

Selected  pre-test  predictions  will  be  discussed  and 
compared  with  post-test  predictions.  The  post- test 
prediction  results  from  the  modification  of  some  of  the 
parameters  in  the  prediction  codes. 

3  EXPERIMENTAL  DATA 

The  test  data  used  here  for  comparison  were  obtained 
in  the  DNW  (Refs.  14,  15).  The  model  is  a  forty  per¬ 
cent  Mach  scaled  model  of  a  hingeless  BO- 105  main 
rotor.  The  blades  are  well  instrumented  and  one  blade 
was  installed  with  124  miniature  pressure  transducers 
and  32  strain  gauges.  The  blade  pressure  and  strain 
gauge  signals  were  recorded  at  a  rate  of  2048/Rev.  An 
average  of  04  revolutions  was  sampled  for  each  rotor 
condition  and  HHC.  The  rotor  performance  data  were 
measured  at  a  rate  of  128  per  revolution  over  a  pe¬ 
riod  of  32  revolutions.  The  wake  geometry  and  vortex 
strength  were  measured  by  the  LLS  and  LDV  tech¬ 
niques  respectively. 

Three  cases  are  selected  for  comparison.  One  is  the 
baseline  case  (Run  140)  without  HHC.  The  second  case 
(Run  138)  produces  lower  noise  with  3/Rev  HHC  (0.8° 
amplitude  and  300°  phase  angle)  and  the  third  one 
produces  (Run  133)  lower  vibration  with  3/Rev  HHC 


(0.8°  amplitude  and  180°  phase  angle).  Other  flight 
parameters  axe: 

Advance  ratio:  /i  =  0.15 
Tip  Mach  number:  Mh  =  0.641 
Shaft  angle:  a,  =  5.3° 

Thrust  coefficient:  Cr  =  0.0044. 

4  RESULTS  AND  DISCUSSIONS 

4.1  Pre-test  Predictions 

Before  the  test  entry  at  DNW,  each  organization  cal¬ 
culated  several  selected  test  cases  with  different  HHC 
settings  and  shaft  angles.  These  results  were  used 
as  a  guide  to  identify  the  BVI  location  for  LLS  flow 
visualization  and  LDV  measurements  as  well  as  a 
sanity  check  for  the  pressure  data.  Although  there 
are  extensive  results  available  for  comparison,  only 
the  predicted  sectional  lift  coefficients  for  the  above 
mentioned  three  cases  are  presented  here.  Fig.  1 
shows  the  comparison  of  the  sectional  lift  coefficients 
at  r/R  =  0.87  for  the  baseline,  low  noise  and  low  vi¬ 
bration  case  respectively.  For  the  baseline  case  (Run 
140),  there  are  eight  BVI  interactions  on  the  advancing 
side  and  five  on  the  retreating  side.  For  the  low  noise 
case  (Run  138),  only  two  dominant  BVI  occurs  on  the 
advancing  side  and  three  on  the  retreating  side.  The 
2/Rev  and  3/Rev  loadings  which  will  affect  the  low 
frequency  component  of  the  radiated  noise  are  best 
predicted  by  ONERA.  Langley’s  calculation  shows  the 
correct  number  of  interactions  on  both  the  advancing 
and  the  retreating  sides  with  higher  magnitude  than 
the  test  data.  It  indicates  that  the  predicted  miss- 
distance  between  the  blade  and  the  vortex  is  probably 
closer  than  the  actual  one.  The  effect  of  blade  dynam¬ 
ics  on  the  blade  BVI  loading  has  been  discussed  in  Ref. 
9. 

4.2  Post-test  Predictions 

After  the  test  was  completed,  each  organization  re¬ 
searched  possible  improvement  to  the  prediction. 
Since  the  dynamic  model  required  improvements  to 
match  the  basic  2/Rev  or  3/Rev  loadings,  AFDD  mod¬ 
ified  the  blade  elastic  properties  at  the  blade  root  sec¬ 
tions  and  used  the  flutter  analysis  in  C  AMR  AD/ J  A 
to  match  the  measured  first  two  flapping  modes,  first 
two  torsion  modes,  and  first  lead-lag  mode.  DLR  at¬ 
tempted  to  improve  their  prescribed  wake  model  by  (1) 
the  measured  vortex  core  radius  from  the  test,  (2)  a 
wake  contraction  of  15  percent,  (3)  adjusting  the  wake 
orientation  with  respect  to  the  blade  tip  path  plane  to 
match  the  measured  vortex  miss-distance,  and  (4)  a 
dual  vortex  pair  consisting  of  a  clockwise  and  a  coun¬ 
terclockwise  rotating  vortex.  NASA  Langley  imple¬ 
mented  the  measured  blade  motion  (flapping,  torsion, 
and  lead-lag)  into  CAMRAD.Mod  1  code.  ONERA 
used  a  measured  vortex  core  size  for  their  calculations. 
The  new  results  are  compared  with  the  test  data. 
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4.3  Tip  Deflection  and  Elastic  Pitch 

The  predicted  tip  deflections  are  plotted  against  the 
measured  results  in  Fig.  2.  Both  DLR  and  ONERA 
results  show  a  reasonable  agreement  with  test  data  for 
all  three  cases  while  AFDD  predictions  do  not  have 
the  2/ Rev  or  3/Rev  shape  even  when  the  modified 
structural  properties  are  used.  It  should  be  noticed 
that  the  difference  of  the  tip  deflection  for  low  noise 
case  at  ^  =  50*^,  where  the  blade  encoimters  the  vor¬ 
tex,  and  rp  =  150®,  where  the  vortex  is  generated,  is 
greater  than  the  one  of  baseline  case.  This  indicates 
that  the  miss-distance  between  the  blade  and  the  vor¬ 
tex  increases  for  the  low  noise  case.  As  for  the  elastic 
pitch  at  =  60®  for  all  three  cases  in  Fig.  2,  most  pre¬ 
dictions  are  over-estimated  as  compared  with  the  test 
data  derived  from  the  measured  strain  gauges.  NASA 
Langley  uses  the  measured  blade  motion  as  an  input 
to  CAMRAD.Mod  1  code,  so  the  results  match  the 
test  data. 

4.4  BVT  Locations  (Tbp  View) 

F^om  the  leading-edge  pressure  measurements  along 
the  rotor  span,  one  may  use  the  valley/peak  at  each 
advancing  and  retreating  side  to  identify  the  measured 
BVI  locations  (Ref.  16).  The  comparison  between  the 
predicted  and  measured  BVI  locations  is  shown  in  Fig. 
3.  For  the  baseline  case  (Rim  140),  the  measured  BVI 
locations  (open  symbol)  occur  in  the  first  and  fourth 
quadrants.  The  measured  BVI  locations  near  ^  =  0® 
seem  to  be  parallel  to  the  flow.  It  is  the  interaction  of 
the  vortex  with  the  rotor  hub  that  causes  the  distortion 
of  the  wake.  For  the  low  noise  case  (Run  138),  the 
measured  BVI  locations  are  clustered  near  %p  =  15®, 
ip  =  90®,  and  ^  =  285®.  As  for  the  low  vibration  case 
(Run  133),  the  measured  BVI  locations  occur  near  rp  = 
60®  and  jp  =  300®.  Only  a  part  of  the  predicted  BVI 
locations  (solid  symbol)  among  eadi  organization  and 
for  all  three  cases  agrees  well  with  the  measured  ones. 
The  difference  in  the  top  view  of  the  BVI  locations 
affects  the  phase  shift  in  the  sectional  loadings.  A 
fifteen  percent  wake  contraction  used  by  DLR  gives  a 
quite  different  BVI  location  near  rp  =  90®  and  rp  = 
270®  as  compared  other  predictions. 

4.5  Wake  Geometry 

The  pre-test  predicted  wake  geometry  and  the  mea¬ 
sured  pressure  data  were  used  to  identify  the  most  im¬ 
portant  blade-vortex  interactions.  However,  the  wake 
geometry  was  measured  before  the  actual  BVI  because 
the  wake  will  be  obscured  by  the  blade.  Selected  vor¬ 
tex  wake  segments  (only  the  tip  vortices  are  presented) 
at  rp  =  35®  and  rp  =  295®  for  all  three  cases  are  com¬ 
pared  with  the  test  data  obtained  from  LLS  technique. 
Fig.  4  shows  both  the  top  and  side  view  of  the  wake 
geometry.  All  the  vertical  coordinates  of  the  wake  ge¬ 
ometry  are  relative  to  the  blade  tip.  No  exact  blade 
position  is  shown  in  Fig.  4,  but  the  blade  is  close  to  a 
horizontal  straight  line  at  z  =  0.  For  the  baseline  case, 


the  miss- distance  prediction  (open  symbol)  of  vortices 
number  5  and  6,  from  each  organization  agrees  well 
with  the  measurements  (solid  symbol)  on  both  the  ad¬ 
vancing  and  retreating  sides,  but  the  top  view  shows 
that  the  interactions  occur  at  different  azimuthal  an¬ 
gles.  The  measured  vortex  segments  were  produced  at 
azimuthal  angles  from  about  rp  =  120®  to  rp  —  150®. 
For  the  low  noise  case,  the  miss-distance  at  ^  =  35® 
is  larger  compared  with  the  one  of  the  baseline  case. 
It  implies  that  the  BVI  loading  is  relatively  weaker. 
Most  predictions  show  an  increase  in  miss-distance. 
For  the  low  vibration  case,  the  blade  produces  nega¬ 
tive  loadings  near  the  outboard  sections.  It  also  pro¬ 
duces  two  trailing  vortices  of  opposite  sign.  One  is 
clockwise  (solid  symbol  with  solid  line)  and  the  other 
is  counterclockwise  (solid  symbol  and  dash  line).  None 
of  the  analyses  predict  this  phenomenon  at  the  present 
time.  However,  a  new  wake  modeling  by  ONERA  that 
deals  with  the  multiple  wake  roll-up  will  be  presented 
in  another  paper. 

4.6  Sectional  Loads 

The  improved  predictions  of  the  sectional  lift  coeffi¬ 
cients  for  all  three  cases  at  r/R  =  0.75, 0.87,  and  0.97 
are  shown  in  Fig.  5.  Predicted  sectional  loadings  by 
AFDD  show  no  significant  improvement  over  the  pre¬ 
test  results.  They  do  not  produce  the  2/Rev  loading 
for  the  baseline  case  and  the  3/Rev  loading  for  the 
other  two  cases.  The  post-test  predictions  by  DLR  do 
show  a  great  improvement  at  the  inboard  sections  for 
all  three  cases.  However,  there  exist  two  loading  spikes 
near  rp  =  90®  and  rp  =  270®  at  r/R  =  0.97.  These  two 
spikes  may  be  caused  by  the  fifteen  percent  wake  con¬ 
traction  that  causes  the  wake  to  nK>ve  inboard  and  pro¬ 
duces  a  large  upwash  outboard.  The  adjustment  of  the 
wake  orientation  also  reduces  the  BVI  peak-to-peak 
loadings  for  some  cases.  The  new  results  predicted  by 
NASA  Langley  are  in  better  agreement  with  test  data 
even  though  some  of  the  2/Rev  and  3/Rev  loadings  are 
over-estimated.  This  shows  that  the  airloads  depend 
strongly  upon  the  blade  motion.  ONERA ’s  results  us¬ 
ing  the  measured  vortex  core  size  show  a  marginal  im¬ 
provement  (some  of  the  BVI  peak-to-peak  values  are 
reduced)  over  the  already  good  pre-test  results. 

6  CONCLUDING  REMARKS 

The  pre-test  results  were  useful  during  the  test  to  lo¬ 
cate  the  important  BVI  events  and  to  reduce  some  of 
the  runs  in  the  test  matrix.  The  HART  test  data  show 
that  BVI  noise  can  be  reduced  with  HHC  (Runl38). 
For  this  case,  the  miss- distance  between  the  blade  and 
the  vortex  wake  was  increased  by  twice  the  value  of  the 
baseline  case  (Run  140).  The  change  in  miss-distance  is 
the  major  factor  in  reducing  the  BVI  noise.  Based  on 
this  fact,  several  alternate  ideas  were  used  to  improve 
the  prediction.  The  wake  geometry  is  improved  by  pre¬ 
scribing  the  motion  of  the  blade  and  this  appears  to 
be  the  most  important  factor  in  obtaining  better  oor- 
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relation  with  the  data.  Other  parameters  such  as  the 
blade  structural  properties,  the  vortex  core  size,  and 
the  wake  orientation  seem  to  be  less  important.  The 
comparison  of  four  different  BVI  prediction  method¬ 
ologies  with  the  HART  test  data  indicates  that  accu¬ 
rate  predictions  require  further  improvements  in  wake 
modeling. 
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Fig,  3  Comparison  of  BVI  locations  (top  view)  with  data  derived  from  measured  leading-edge  pressure. 
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Fig.  4  Comparison  of  measured  and  predicted  vortex  segments  at  rp  =  3S  degrees  and  ^  =  295  degrees. 
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Fig.  4  (continued)  Comparison  of  measured  and  predicted  vortex  segments  aX  xp  =  35  degrees  and  xp  =  295  degrees. 
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Fig.  4  (continued)  Comparison  of  measured  and  predicted  vortex  segments  at  V’  =  35  degrees  and  V"  =  295  degrees. 


Fig.  4  (continued)  Comparison  of  measured  and  predicted  vortex  segments  at  ^  =  35  degrees  and  ^  =  295  degrees. 


